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ABSTRACT 

Using the Australia Telescope Compact Array we have detected CO (1-0) 
and CO (5-4) from TN J0924-2201 at z = 5.2, the most distant radio galaxy 
known to date. This is the second highest redshift detection of CO published 
so far. The CO (1-0) is 250 — 400km s -1 wide with a peak flux density of 
520 ±115 /iJy beam -1 whilst the CO (5-4) line emission is 200 — 300km s -1 wide 
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with a peak flux density of 7.8 ±2.7 mJy beam 1 . Both transitions are spatially 
unresolved but there is marginal evidence for spatial offsets between the CO and 
the host galaxy; the CO (1-0) is located 28 ± 11 kpc (4"5 ± V!7) north of the 
radio galaxy whilst the CO (5-4) is located 18 ± 8 kpc (2"8 ± 1"2) south of the 
radio galaxy. Higher spatial resolution observations are required to determine 
the reality of these offsets. Our result is the second detection of CO in a high 
redshift galaxy without pre-selection based on a massive dust content. 

Subject headings: early universe — galaxies: formation — galaxies: high redshift 
— galaxies: individual (TN J0924-2201) — radio lines: galaxies 

1. Introduction 

Right out to the highest redshifts so far identified, radio galaxies are the most massive 
systems at each epoch. They reside in regions of over-density (Le Fevre et al. 1996; Carilli 
et al. 1997; Pascarelle et al. 1998; Pentericci et al. 2000), corroborating evolutionary models 
which trace high redshift radio galaxies (z > 2; HzRGs) through to nearby central dominant 
elliptical galaxies in cluster- type environments (e.g. Lilly & Longair 1984). Unlike their 
nearby counterparts, however, HzRGs are not massive ellipticals with old stellar populations. 
Instead, their rest-frame optical/UV morphologies are composed of multiple blobs of clumpy 
emission often aligned along the axis of relativistic jets (van Breugel et al. 1998; Pentericci 
et al. 1999). Rest frame far infra-red observations are strongly suggestive of large-scale, 
massive, dust-obscured star formation (Stevens et al. 2003). Furthermore, HzRGS are often 
enshrouded by massive gas reservoirs which extend over tens of kpc, traced both by Lyman- 
a halos (McCarthy et al. 1990; Villar-Martm et al. 2003) and by carbon monoxide (CO ) 
line emission (Papadopoulos et al. 2000; De Breuck et al. 2003). Each of these properties is 
consistent with the scenario where the galaxies are being observed in their formative stages. 
It is essential, therefore, that the characteristics of the HzRG population are investigated in 
detail in order to further our understanding of galaxy, and cluster, formation and evolution. 
In this Letter we report the detection of molecular gas at z — 5.2 in TN J0924-2201, the 
highest redshift radio galaxy and the second highest redshift CO detection to date, after the 
z = 6.42 quasar SDSS J1148+5251 (Fan et al. 2003; Walter et al. 2003). Where necessary 
we use a flat A-dominated cosmology with H = 71 km s _1 Mpc -1 and f2 A = 0.73 (Spergel 
et al. 2003). In this cosmology, 1" = 6.3 kpc at z — 5.2. 
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2. TN J0924-2201: selection & characteristics 

TN J0924-2201 was initially discovered (De Breuck et al. 2000) from radio surveys using 
the well established technique for finding HzRGs of selecting ultra-steep spectrum sources 
with faint K-band counterparts (e.g. Chambers et al. 1996; Blundell et al. 1998). Follow-up 
optical spectroscopy detected faint Lyman-a emission from the host galaxy at a redshift 
of 5.1989 ± 0.0006 (van Breugel et al. 1999; Venemans et al. 2004). Based on its compact 
(1'.'3) radio morphology, under-luminous Lyman-a emission compared with other HzRGs and 
rest frame 350 nm emission composed of several clumpy optical features, it was suggested 
that TN J0924-2201 is a young, primeval galaxy actively forming the bulk of its stars (van 
Breugel et al. 1999). We initially chose TN J0924-2201 as a candidate for CO observations 
because, as the most distant known member of a sample of HzRGs, it is presumably one 
of the youngest (Blundell & Rawlings 1999). The existence of molecular gas would provide 
independent evidence for ongoing massive star formation, expected if we are probing its 
formation epoch. The southern declination (5 = —22°) of this source makes observations 
from millimetre telescopes in the northern hemisphere difficult, and its redshift moves all 
CO transitions out of any Very Large Array (VLA) band. We used the new wide band- 
width 1.2 cm and 3 mm receivers on the Australia Telescope Compact Array (ATCA) to 
search for CO (1-0) and CO (5-4) with z/ rest = 115.2712 GHz and 576.2677 GHz, respectively. 



3. ATCA Observations 

Observations were carried out with the ATCA during 2004 August and September. The 
ATCA is an array of six 22-m antennas with baselines ranging from 30 m to 6 km and capable 
of East- West, North-South and Hybrid array configurations. We utilised the maximum 
available bandwidth of the current correlator with 128 MHz/64 channel configurations for 
two adjacent IF bands. The spectral resolution of this configuration is 4.4 MHz. For the CO 
(1-0) observations, we used the hybrid H 168 array 1 , tuning the IFs to 18.561 and 18.631 GHz 
respectively, allowing an overlap region (after bandpass correction) of 20 MHz. This setup 
benefited from a simultaneous velocity coverage of 2500 km s _1 and resolution of 72 km s _1 . 
System temperatures ranged from 30-50 K and phase stability was within ±20° over the 
course of the observations. Phase, amplitude, bandpass and astrometric calibration were 
acquired by short observations of PKS B0919-260 every 12-25 minutes depending on the 



1 Thc size of the naturally weighted restoring beam in the H 168 array at 18.6 GHz is 14'.'5 x lO'.'l, or 
90 x 63 kpc. 
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weather. Absolute flux density calibration was ensured by snapshot observations of Mars. 
We estimate the calibration to be accurate to better than ±20%. 

The observational setup was repeated in the H 75 array 2 with the IFs tuned to 92.853 
and 92.921 GHz respectively for the CO (5-4) transition. System temperatures (measured 
once per hour using a 300 K mechanical paddle moved over the feedhorn) ranged from 
200-350 K and phase stability was mostly within ±50° over the course of the observations. 
Phase and amplitude calibration were based on observations of PKS B0834-201 every 12 
minutes and the bandpass was determined by 10 minute observations of either PKS B1253- 
055 or PKS B0537-441. The rms error in the pointing ranged from 3 — 7" depending on the 
antenna (the primary beam at 93 GHz is 31"). We periodically observed the 0.5 Jy quasar 
PKS B0925-203, less than 2° from TN J0924-2201, to verify the calibration and astrometry. 
We estimate the flux scale to be accurate to better than ±30%. Data reduction was carried 
out in accordance with standard ATCA calibration procedures for millimetre observing 3 . 
Overall, 40 hours and 35 hours of observing time were employed at 18.6 GHz and 92.9 GHz 
respectively to carry out this detection experiment. Table 1 gives a summary of our results. 



4. Results 
4.1. CO (1-0) 

The CO (1-0) line is detected in at least five independent channels covering ~ 400 km s -1 and 
shown in Figure 1. We subtracted a 0.71 mJy continuum in the UV plane using a first order 
spectral baseline fit. The continuum image is shown in greyscale in Figure 2. The velocity- 
integrated, continuum-subtracted, flux density is 0.087 ± 0.017 Jy beam -1 km s" 1 and the 
corresponding intrinsic CO line luminosity is L' co = 1.2±0.26 x 10 11 K km s -1 pc 2 . There are 
two kinematic features apparent in the spectrum: the brighter, dominant, emission centred 
at v = 150 km s -1 (z = 5.202 ± 0.001) and a weaker feature lying ~ 200 km s" 1 redward. 
Deeper observations are required to determine if the weaker feature is real. As commonly 
observed, the CO lies redward of the Lyman-a derived redshift (see Table 1 in Klamer 
et al. 2004). There is no conclusive evidence for spatially extended gas, but the velocity- 
integrated emission (Figure 2) appears to be offset to the north of the radio galaxy by 
4'.'5 ± 1'.'7 (28 ±11 kpc). At this stage we do not consider this offset to be significant, since 



2 The size of the naturally weighted restoring beam in the H 75 array at 92.8 GHz is 5'.'9 x 4'.'7, or 
37 x 30 kpc. 

3 ht t p : / / www . at nf . csiro.au/ computing/miriad 
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it is much less than a beam size away from the radio galaxy position. For comparison, the 
position angle of the projected radio axis is 29° (see Figure 2). Deeper observations with 
higher spatial resolution are required to better constrain the spatial extent of this gas. 

4.2. CO (5-4) 

At 93 GHz our instantaneous velocity coverage was only 600 km s _1 with a spectral 
resolution of 15 km s -1 . The resultant spectrum is shown in Figure 1 and the emission 
profile integrated from to 300 km s" 1 is shown in Figure 2. Although the zero level in the 
profile is well determined from offset positions in the image, our bandwidth is only marginally 
wider than the line width, so we have to consider the possibility of a spectral baseline offset 
due to continuum emission from the radio galaxy. However, the flux density we observe is 
more than 100 times that expected from the continuum. This is illustrated by the radio 
to sub-millimetre spectral energy distribution shown in Figure 3. Extrapolating the non- 
thermal synchrotron spectrum with a spectral index of a = —1.7, we expect a contribution 
of < 0.04 mJy whilst extrapolating from the 352 GHz SCUBA limit assuming a modified 
blackbody spectrum with j3 = 1.5 and T=50 K, the expected contribution due to thermal 
dust radiation is < 0.01 mJy. The most likely source of the 93 GHz emission is therefore 
the redshifted CO (5-4) transition at z — 5.202. The velocity-integrated flux density is 
1.19 ± 0.27 Jy beam -1 km s _1 , corresponding to L f co = 4.3 ± 1.0 x 10 10 K km s" 1 pc 2 . 
Again, the emission appears offset from the radio galaxy position, the peak lying 278 ± 172 
(18 ±8 kpc) south. Again, the positional offset is less than a beam size. Deeper observations 
are required to constrain the spatial extent of the molecular gas, including whether there are 
physically distinct gas emitting regions. 

5. Physical properties of the molecular gas 

5.1. Mass Estimates 

CO observations are translated into total H 2 mass estimates, assuming that the observed 
CO traces molecular clouds on large scales and that the molecular gas mass dominates the 
total dynamical mass of the galaxy. These mass estimates are uncertain at best, but to 
facilitate comparison with observations from other HzRGs, we assume M(H 2 ) / L' co = a = 0.8 
M (K km s _1 pc 2 ) -1 as measured for nearby ultraluminous infra-red galaxies (Downes & 
Solomon 1998). Since molecular clouds provide both the fuel and the site for the next 
generation of star formation, our result implies that there is ~ 10 11 M of available gas 
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from which to manufacture the next stellar generation. This estimate is consistent with 
M(H 2 ) estimates of other HzRGs like 4C 60.07 (Greve et al. 2004; Papadopoulos et al. 2000) 
and B3 J2330+3927 (De Breuck et al. 2003). There is no direct evidence to indicate that 
TN J0924-2201 is gravitationally lensed. However, if lensing is important then L' co , and 
hence M(H 2 ), are upper limits only. 

5.2. Temperature and Density from LVG modelling 

The ratio of the peak flux densities would be equal to the square of the ratio of the 
upper transition levels (1/25) in the event that both lines are optically thick and thermally 
excited. Under the assumption of optically thick gas, the ratio we measure is 1/(15 ± 9), 
which suggests the CO (5-4) transition is sub-thermal; however, the errors are large. We 
have used a standard single-component large velocity gradient (LVG) code, with a typical 
CO-to-H 2 abundance X[CO]/(dv/dr)=8 x 10" 5 (km s" 1 ) _1 pc, to further constrain the 
physical properties of the molecular gas. A good match to the observations is given by 
log(n(H 2 )) = 3.3 (where n(H 2 ) is measured in units of cm" 3 ), T kin = 50 K and r> 2.0 kpc, 
but other density/temperature combinations are also consistent with the observed fluxes (see 
Figure 4). Observations of CO transitions higher than CO (5-4) are required to constrain the 
excitation conditions of the gas. If the observed offsets between the CO (1-0) and CO (5-4) 
transitions reflect physically distinct CO emitting regions, then the excitation conditions 
presented here become invalid. 

5.3. Where is the dust? 

At present, the most efficient technique for detecting high redshift molecular gas is to 
pre-select the dustiest objects, as determined from their large rest frame far-infra red lumi- 
nosities using (sub) millimetre sensitive instruments like SCUBA on the JCMT or MAMBO 
on the IRAM 30m telescope. This also means that current samples are biased towards ob- 
jects with the largest dust-to-gas ratios (Lfir/^co ~ 300 L K" 1 km" 1 s pc" 2 ; Carilli 
et al. 2004). For comparison, the median ratio is 160 for nearby ultra-luminous infra red 
galaxies with LpiR > 10 12 L (Solomon et al. 1997) and about 50 for nearby luminous infra 
red galaxies with 10 n L Q < Lfir < 10 12 L o (Gao & Solomon 2004). The selection bias at 
high redshift is further compounded by the fact that initial CO detections are usually made 
in J > 3 order lines, which prejudices the sample against regions of low density molecular 
gas where the higher order transitions are sub-thermally excited. 
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Five of the six published CO detections of HzRGs constitute the most FIR-luminous 
z > 1 radio galaxies recently surveyed with SCUBA (Reuland et al. 2004). In this sample, 
approximately 50% of the 24 observed galaxies were not detected; TN J0924-2201 featured 
amongst these non-detections with L 850fim < 1 x 10 24 W Hz -1 , five times underluminous 
compared to the CO detected HzRGs in the sample. Therefore, TN J0924-2201 will have a 
dust-to-gas ratio at least five times smaller than the current sample of high redshift galaxies, 
similar to the nearby luminous infra-red galaxy population. The only other CO emitter 
at high redshift not selected on the basis of FIR emission — that we are aware of — is 
53W002 at z = 2.394 (Yamada et al. 1995; Scoville et al. 1997; Alloin et al. 2000). This 
source has a 3a upper limit of 4.3 mJy at 850/zm (Archibald et al. 2001), corresponding to 
^850 M m < 7x 10 24 W Hz" 1 . An a posteriori detection of the dust content in 53W002 (1.7±0.4 
mJy at 1.3 mm; Alloin et al. 2000) suggests that a similar observation of TN J0924-2201 
may also reveal its dust content. Thus 53W002 and TN J0924-2201 may be representative 
of a population of high redshift sources which remain undetected due to current selection 
techniques. 



6. The co-evolution of a central Dominant and a protocluster 

At z — 5.2 the universe is only about a billion years old — time enough, however, to 
build a supermassive black hole, trigger powerful radio jets and produce a stellar population 
capable of enriching the interstellar medium with ~ 10 7 M of CO (assuming an abundance 
of CO relative to H 2 of Zco = 10~ 4 ). This rapid enrichment could be facilitated in part by 
the powerful radio jets themselves propagating into overdensities in the primordial gas, and 
shock-inducing massive star formation (Klamer et al. 2004; Bicknell et al. 2000). 
Venemans et al. (2004) recently postulated that TN J0924-2201 resides within a protocluster 
environment after discovering six Lyman-a emitters in a 6'8 x6f8 field (2.6 x2.6 Mpc) around 
the radio galaxy — the typical scale size for nearby rich clusters (Abell 1958). Our CO 
observations substantiate a five-year old hypothesis that TN J0924-2201 is a young forming 
galaxy, potentially the primeval central dominant galaxy of a rich cluster. Thus it seems 
that cDs may form with their host environments from the earliest times in the Universe. 

Big thanks go to everyone who made the 12 and 3 mm ATCA upgrades possible. Thanks 
also to Chris Carilli for the VLA observations which we have incorporated into Figures 2 
and 3 and Christian Henkel who originally wrote the LVG code. UK acknowledges receipt 
of an Australian Postgraduate Award and CSIRO Postgraduate Scholarship. The Australia 
Telescope Compact Array is part of the Australia Telescope which is funded by the Com- 
monwealth of Australia for operation as a National Facility managed by CSIRO. 
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Fig. 1. — Spectra of the continuum-subtracted CO (1-0) (shaded region) and CO (5-4) (dashed 
line) transitions, through the peak of their respective velocity-integrated emission shown in Figure 2. 
The CO (1-0) spectrum is unbinned (32 km s _1 bins; spectral resolution 72 km s _1 ) whilst the 
CO (5-4) spectral channels have been re-binned into 30 km s _1 bins (original spectral resolution 
15 km s _1 ). The velocity axis is defined with respect to the Lyman-a derived redshift of z = 5.1989 
and the bars in the top left and right corners show the rms noise level in the CO (1-0) and CO 
(5-4) cubes respectively. 
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Fig. 2. — Greyscale 18.6 GHz continuum image with CO (1-0) and CO (5-4) velocity-integrated 
emission overlaid in red and cyan respectively. The continuum flux density is 0.71 ± 0.03 mJy, the 
CO (1-0) contour levels are ±2,3,4,5 x a where a = 0.017 Jy beam _1 km s _1 and the CO (5-4) 
contour levels are ±2.5,3,3.5,4 x a where a = 0.27 Jy beam" 1 !™ s" 1 . The sizes of the 18 GHz 
(open ellipse) and 93 GHz (filled ellipse) restoring beams are shown in the bottom left corner, the 
cross shows the optical position of the radio galaxy. The radio lobes, illustrated with filled circles, 
extend l'.'3 along a position angle of 29°. 
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Fig. 3. — Radio to sub-millimetre spectral energy distribution of TN J0924-2201. The stars are 
radio observations from the Texas, NVSS and PMN radio sky-surveys along with pointed VLA 
observations and the ATCA 18.6 GHz continuum point from this work. The solid line is a power 
law with spectral index a = — 1.7. The 352 GHz upper limit is from Reuland et al. (2004) and the 
dashed line shows the Rayleigh-Jeans portion of a modified blackbody spectrum, normalised to the 
352 GHz (850 fim) limit, with T=50 K and = 1.5. The upper limit at 93 GHz is from this paper 
using the CO (5-4) flux density averaged over the range 0-300 km s _1 . 
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Fig. 4. — CO ladder with LVG models overlaid. The best match to the observations is 
log(n(H2)) = 3.3, Tki n = 50 K and r> 2.0 kpc (solid line). Also in agreement with the obser- 
vations are log(n(H2)) = 4.0, Tk in = 30 K and r> 2.9 kpc (dashed line) and log(n(H2)) = 2.7, 
Tkin = 150 K and r> 2.2 kpc. (dotted line). The condition for optically thick, thermally ex- 
cited gas is the solid, parabolic fit. The error bars are la noise estimates and do not include the 
uncertainty in the absolute flux scale (§3) 
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Table I. Observed physical parameters of TN J0924-2201 





18.6 GHz 


CO (1-0) 


CO (5-4) 




continuum 


line 


line 


Peak (mjy beam 1 ) 


0.71 ±0.03 


0.52 ±0.12 


7.8 ±2.7 


Width (km s- 1 ) 




250 - 400 a 


200 - 300 


I C o (Jy km s _1 ) 




0.087 ±0.017 a 


1.19 ±0.27 


ARA J2000 (") b 


-0.7 ± 0.7 


1.5 ±2.1 


-1.0 ±1.3 


ADEC J20 oo(") b 


0.3 ±0.7 


4.5 ±1.7 


-2.8 ±1.2 



Note. — "The weaker feature in the CO (1-0) spectrum has 
not been included in this measurement because within the 300-400 
km s _1 velocity range, there are several other similar features in 
the image which decrease its significance. If this feature were in- 
cluded the velocity-integrated flux density of the CO (1-0) emis- 
sion would be 0.12 ± 0.03 Jy km s _1 .''These positions are offsets 
from high resolution VLA observations: RA J20 oo = 09 h 24 m 19.92 s , 
DEC/2000 = -22°01'41.5" (van Breugel et al. 1999). Positive offsets 
are North and East of the radio galaxy position. 



